General Disclaimer 


One or more of the Following Statements may affect this Document 


• This document has been reproduced from the best copy furnished by the 
organizational source. It is being released in the interest of making available as 
much information as possible. 


• This document may contain data, which exceeds the sheet parameters. It was 
furnished in this condition by the organizational source and is the best copy 
available. 


• This document may contain tone-on-tone or color graphs, charts and/or pictures, 
which have been reproduced in black and white. 


• This document is paginated as submitted by the original source. 


• Portions of this document are not fully legible due to the historical nature of some 
of the material. However, it is the best reproduction available from the original 
submission. 


Produced by the NASA Center for Aerospace Information (CASI) 




Semi-Annual Status Report 
to the 

NATIONAL AERONAUTICS AND SPACE ADMINISTRATION 

on 

NASA Grant NSG-3048 
ALTERNATIVES FOR JET ENGINE CONTROL* 
October 1, 1982 - March 31, 1983 


*This work has been carried out under the direction of 


Dr. Michael K. Sain 
Department of Electrical Engineering 
University of Notre Dame 
Notre Dame, Indiana 
USA 46556 


ABSTRACT 


This report deals with progress made on the Grant NSG-3048 during the 
six month period beginning October 1, 1982 and ending March 31, 1983. The 
NASA Technical Officer for this period was Dr. Kurt Seldner of Lewis Re- 
search Center, The director of the research at the University of Notre 
Dame was Dr. Michael K. Sain, who has been assisted by Mr, Daniel Bugajskl 
and Mr. Joseph O'Sullivan, graduate research assistants in the Department 
of Electrical Engineering. Mr. Thomas A. Klingler has received the degree 
of Master of Science based upon his studies under this Grant, for his 
thesis entitled "Nonlinear Modeling of a Turbofan Jet Engine: A Tensor 

Method Application". Mr. Klingler is presently with the Pontiac Division 
of General Motors Corporation, where he is a principal in the moderniza- 
tion of the Division’s automobile engine simulation and test facilities. 

Researches during this period have stressed tensor model order re- 
duction, recursive tensor model identification, input design for tensor 
model identification, software development for nonlinear feedback control 
laws based upon tensors, and development of the CATNAP software package 
for tensor modeling, identification, and simulation. This report is 
focused upon the last of these. The other projects are scheduled for 
future reports. 
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SECTION I 


INTRODUCTION 

Models have always been an important aspect of applications engine- 
ering in the area of multivariable control system design. See for ex- 
ample [13. Practical industrial examples of the use of models are pro- 
vided by gas turbine engine control systems, which commonly use models 
to generate control and response trajectories for various power demands. 
These models, when scheduled over the operating envelope, can reduce the 
compensation normally required of the controller, and thus provide the 
feedback loop with an opportunity to achieve better accuracy in the pre- 
sence of noise and parametric uncertainties. 

The scheduling of feedforward models and feedback compensation typ- 
ically produces nonlinearities, even if the local models are linear. 
Accordingly, there is a basic interest in fundamental approaches which 
incorporate nonlinearity at the outset. Such approaches should (1) re- 
duce to the earlier schemes for variables with very small excursions, 

(2) be amenable to scheduling, and (3) offer opportunities for determin- 
ation of parameters from simulation data. 

One such approach uses the notions of power series and algebraic 
tensors [2] to generate a class of nonlinear models [3]. The impor- 
tant feature of the algebraic tensor is that it provides an organized 
way of describing the power series expansion formula, lending itself 
with relative ease to programming on a digital computer. 

This report summarizes an investigation focused on the use of the 
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above tecUtlque for the Identification and testing of the algebraic non- 
linear model structure which utilizes the tensor approach [4]. The 
basic strategy involves initial condition and control signal design in 
order that the nonlinear mode] not only outperforms a standard linear 
model locally about an expansion point, but also establishes a larger 
region of model validation. Emphasis throughout the study is placed on 
the application of these ideas in nonlinear model generation using a 
turbofan jet engine simulation. 

The remainder of the report is outlined as follows. In Section II 
the background material necessary for an understanding of the modeling 
methodology is presented. Section HI describes the modeling software 
which has evolved into a computer-aided design package capable of handl- 
. ing the overall process of model formulation, identification, and vali- 
dation for any particular nonlinear modeling pv>-. .'m developed within 
the tensor framework. Next, Section IV presents an overview of the QCSE 
engine along with its corresponding digital simulation. QCSEE is the 
name of the turbofan jet engine used in the modeling exercises of this 
investigation. A detailed identification with full validation studies, 
as well as a summary of the modeling technique as applied to QCSEE, is 
included in Section V, Finally, Section VI offers a few pertinent 
closing remarks. 


SECTION TI 


NONLINEAR MODEL STRUCTURE 

The purpose of this section is to relate some of the background ma- 
terial which is fundamental to the work described in the succeeding sections . 
The first portion summarizes important ideas developed in the previous 
work [3 ] which serves as a basis for this study. Since a rigorous de- 
velopment of these topics is not the purpose of this report, the presen- 
tation here is limited. However, the discussion should supply an adequate 
foundation for the material to follow. Finally, the section closes with 
the introduction of a technique used to normalize the bulk data required 
in the identification phase, , 

2.1 ALGEBRAIC TENSORS IN MODELING 

Current research indicates that the use of tensor algebra £ 2 J , in 
modeling nonlinear dynamical systems, is quite advantageous. Furthermore, 
these studies maintain that the use of algebraic tensors gives a unified 
view of the nonlinear modeling problem. The Important feature of the 
algebraic tensor is that it provides an organized way of describing the 
power series expansion formula. In addition, its use allows for the im- 
plementation of linear parameter identification techniques. 

Due to the brevity of the following treatment, the reader is advised 
to consult the works [3,5,6] for detailed expositions of the topics 
discussed. With this in mind, consider the following nonlinear ordinary 
differential equation 

x - f(x,u) . (2-1) 
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Here x e X is a n-tuple of states and u e U is a m-tuple of inputs, 
where X and U are real vector spaces. The function f : X MI+X, 
under certain assumptions, may be expressed in a power series expansion 
of two variables. Moreover, if we employ the universal bilinear tensor 
product function » : X x U + X®U [7], as a powerful tool to organize 
the higher order mixed partial derivatives from the expansion, then the 
right side of (2-1) can be expressed in the form 


f(x,u) « J T L. . (x © x ® . . . ® x ® u ® u ® ... ® u). (2-2) 

i-0 j»0 1J ' v— 1 ' v- ' 

i times j' times 

Therefore, a nonlinear model of (2-1) can be obtained by a truncation of 
(2-2). Substituting (2-2) into (2-1), and rewriting the result in a dif- 
ferent form, yields the model structure 

* A/ ^ 

* " EL 10 L 01 L 20 L 11 L 02 "30 


) X* 


(2-3) 


where x^ is a vector stacked with ordered tensor product, terms given 


by the expression 


T 

X T ■ (x, u, x ® x, x ® u, u ® u, x ® x ® x, . . .) . (2-4) 

It should be noted that contains a number of redundant terms due to 

the commutativity of scalar multiplication in the field; however, through 
use of the symmetric tensor algebra [ 7 ] these redundancies are elimi- 
nated and a reduced-size version of (2-3) is constructed. 

<v 

To identify the parameters, (2-1) is perturbed from (x o ,u q ) 

and excited with cosinusoidal inputs. The state solutions and control 
Inputs are sampled at h selected time points and are loaded into the 
pxh matrix X^. The sampled values of x and u comprise the first 
n+m rows of X,p, whereas the remaining p-(n+m) rows contain monomial 
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terms which are products of various combinations of the first n+ra rows. 
These monomial terms are generated via an efficient ordering algorithm [ 3 ] 
which minimizes the required number of multiplications by re-using any pro- 
ducts previously computed instead of repeating the operations. The nxh 
matrix X is formed by loading derivative estimates for x at the h 
points in time. As an illustration, a truncation approximation which re- 
tains tensor product terms through second degree is given by the equation 

X - rt„ x_ . 


10 01 “20 11 02 - 


(2-5) 


Finally, to arrive at the desired operators, a singular value 

decomposition is employed to solve the minimal least-squares problem. 


2.2 DATA NORMALIZATION 


Recall from the previous section that sampled data, in the form of 
derivative estimates, state solutions, and control Inputs, is used to 
construct a least-squares minimization. Frequently though, when physical 
,'ystems are being considered, this raw data spans a wide range of units 
and therefore taxes the abilities of the coefficient matrix software* 
Consequently the need for a practical, yet easy to implement, normaliza- 
tion technique is essential. The rest of this section is devoted to out- 
lining one such scheme. 

To begin, let 

L - '*10 %1 *20 *11 *02 *30 - 1 * 

then by substitution equation (2-3) becomes 

X - LX t . (2-7) 

• • 

Now, let Z * VX, where V is a normalization matrix given by 
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Similarly let Zy ■ WX^. In this situation the normalization matrix W 
takes on a rather complicated form due to the nature of X^,, but through 
use of the KrOnecker product [ 8 ] of two matrices, W can be defined as 
the partitioned matrix 



pxp 
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Applying this transformation to (2-7) yields the normalised matrix equa- 
tion 

Z - L Z T (2-3) 

i 

which can than be solved for L. The data in Z and Z ^ now takes on 
real values between -1.0 and 1.0. 

Work backwards now to reconstruct the original system. Equation 

# 

(2-8) can be written as VX - LVBCj, and manipulated to equal (2-7) when 
L - V _1 LW. 

This approach will prove to be quite valuable when the turbofan 
jet engine problem is tackled in Section V. 

In closing, this section has addressed certain issues of theoreti- 
cal interest which serve as a prerequisite to the material contained in 
if'U? jnges that follow. The remainder of the report centers around the 
application of these ideas. 


SECTION III 


CATNAP 

The intent of this auction is to present a detailed discussion of the 
£omputer-Aided Tensor Nonlinear Modeling Applications Package (CATNAP) cur- 
rently used as a development tool in the formulation, identification, and 
validation of nonlinear models of the type mentioned in Section II. 

3.1 OVERVIEW 

The structure of CATNAP is based upon ideas from distributed proces- 
sing and local networking [ 9 ] in which computations are spread over mul- 
tiple machines. The package [10], developed using the extensive capa- 
bilities of the IBM and DEC Command Procedure Languages with the strengths 
of FORTRAN and SPEAKEASY, is tailored to utilize existing computer hard- 
ware effectively, as well as to provide the fastest possible user turn- 
around time. Development of a package such as CATNAP is an absolute 
necessity to apply the ideas of {3], 

Specifically, CATNAP is divided into two segments, each of which is 
implemented on an independent computing system. These segments are en- 
titled GENERATE and SIMULATE. GENERATE is implemented on a Time Sharing 
Option [11] of the IBM 370-168 mainframe computer system, where the mem- 
ory dependent and highly computational routines of the package can bene- 
fit from use of the virtual memory and floating point hardware. Once a 
structured nonlinear model is identified, it is then transfered to the DEC 
PDP 11/60, where, using SIMULATE, it can be analyzed and compared to a 
linear model as well as the true system. In this manner, the user has 
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at his disposal both the graphics capabilities of a video terminal and an 
electrostatic printer/plotter for the immediate display of data and com- 
parative trajectories. Furthermore, both of these segments are highly in- 
teractive and contain straightforward input prompts as well as informative 
error messages. 

3.2 GENERATE 

The GENERATE segment of CATNAP is primarily made up of three routines 
governed by a higher level supervisor. This supervisory level is written 
using the command procedure language CLIST [12] and performs the following 
main functions in sequence: 

(1) prompts the user for the name of the desired loader rou- 
tine to be executed; 

(2) passes control to the loader routine defined in (1) ; 

(3) passes control to IDENTIFY; and, 

(4) upon user request, passes control to TRANSFER. 

In addition to these main functions, certain maintenance roles such as 
file creation, allocation and deletion are handled by this supervisor. 
Figure 3.1 contains a block diagram depicting the structure of GENERATE. 

Associated with each nonlinear system to be modeled, there exists a 
loader routine which performs the model formulation task. These routines 
are stored in a library and are typically written in double precision 
FORTRAN. The purpose of any loader routine is to excite the given non- 
linear system via initial condition and control input perturbations and 
to sample the states, inputs and derivative estimates over h selected 
points in time. The system is then represented by the matrix equation 
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X " fL 10 L 01 L 20 L 11 L 02 L 30 X T * 

I. 1 L 

nxh nxp pxh 

The first n-Hn rows of the matrix are formed from the sampled values 

of x and u while the remaining p-(n+m) rows are formed from combi- 
natorial products of those first n+m rows. X is formed by loading 
derivative estimates for x^, ^i****^ at the h time points. The 
number p depends on n and m, and the degree of the truncation approx 
imation. All this data is then stored in TEMPFILE for later use. 

The ordering algorithm which is used to construct X^ utilizes the 
important characteristic that higher degree tensor terms are "built up" 
from lower degree terms. Incorporating this idea into the methodology 
reduces the computational aspects of the scheme significantly. In fact, 
this improved algorithm operates nearly 3.5 times faster for second de- 
gree approximations than the previously used method. Even a greater im- 
provement is expected for higher degree approximations. 

Using the above formulation technique, CATNAP can accommodate any 
particular nonlinear modeling problem since the problem specifics are 
transparent to the remainder of the package. The only requirement is 
that TEMPFILE contains the appropriate data, 

After the completion of any chosen loader routine, the program 
IDENTIFY is executed. IDENTIFY reads the interim data from TEMPFILE 
and forms a least squares minimization problem which is solved for the 

r\j r\, 

partitioned matrix containing the desired parameters. These 

parameters are recorded at the terminal as well as entered into the 
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MODEL data file. It should be noted here that IDENTIFY is written in the 
high level language SPEAKEASY [13], which is based on the concepts of ar- 
rays and matrices and processes these as entities. This results in the 
elimination of the many loops necessary in other programming languages. 

The main reason for employing SPEAKEASY here is that the highly efficient 
routine SIMEQUAT can be easily used to solve the least squares problem via 
singular value decomposition, thus reducing the apparent complexity of the 
problem to a minimum. 

Once the identified data is entered into the MODEL data file, the 
user may elect to execute the program TRANSFER. Upon a yes response to 
a supervisory prompt, the program is submitted batch to the IBM 370-168. 
TRANSFER is merely a Job Control Language (JCL) deck which sends a copy 
of the file MODEL, containing the parameters, to the DEC PDP 11/60 

computing system by the way of a Remote Job Entry port, and stores it in 
the nonlinear model library. An excellent account of JCL can be found ^n 

[14], 

3.3 SIMULATE 

Shifting our concern away from the discussion of GENERATE, we now 
focus our attention on the SIMULATE segment of CATNAP. Basically, two 
routines plus a supervisory comprise the structure of SIMULATE. Figure 
3.2 offers an illustration of this structure to supplement the following 
presentation. Written in the form of an Indirect Command File [15] this 
supervisor allows the user to: 

(1) create new simulator routines; 
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(2) execute existing simulator routines; and, 

(3) execute VF.RSATEC which produces hardcopy plots. 

As earlier, this supervisor performs a number of file maintenance duties 
In addition to the above functions. 

Fundamentally, the heart of this segment is contained within any 
given simulator routine. Typically a simulator routine, Written in FORTRAN, 
exists for each nonlinear modeling problem studied; however, only one sub- 
routine in that program is altered among versions, and that is the appli- 
cation subroutine TP.UES. The remainder of the program stays unchanged. 

TRUES contains the true system representation of the nonlinear system 
being modeled, and is used extensively in comparison studies. Because of 
the number of TRUES subroutines that exist, a library has been created to 
store the various simulator routines. 

The execution of a particular simulator routine can be divided into 
three steps; (1) problem configuration; (2) systems integration; and (3) 
solution display. The first of these steps requires the user to decide 
which of the available systems, true solution, linear model and/or non- 
linear model, should be included in the session configuration. When a 
model is chosen, the user is asked to enter the name of the desired model. 
That model is then read into the simulator from the appropriate library. 

The linear models which are used in CATNAP are generally identified by 
standard .techniques and are available for use in comparison studies. 

Next, the user is prompted for various integration parameters such 
as stepsize and upper time limit as well as initial conditions, input 
amplitudes and frequencies. The configured systems are then integrated 
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and the solutions are sampled at 100 points, evenly spaced In time. Fin- 
ally, to assist In the data analysis, a number of options are available 
to the user. They Include: 

(1) printing the solutions on the Versatec; 

(2) displaying the comparative trajectories on the Tek- 
tronix graphics terminal; 

(3) writing the trajectory solutions to SPOOL for hard- 
copy plotting at a later time; and, 

(4) solving the configuration for another set of initial 
conditions and control inputs. 

The use of these options provides a powerful yet flexible capability for 
the study of model performance and validity. Furthermore, when all three 
systems are includedin the configuration, an additional error criterion 
is generated and used in testing the performance of the nonlinear model 
versus that of the linear model (See [3]). 

N 

Let denote the absolute maximum error in the nonlinear model 

solution, as compared to the true simulation solution, over the time 

range of simulation for the ith state variable. Similarly, define 

N L 

for the linear model error. Then is the difference - e^. Thus, 

if is negative, the nonlinear model has exhibited a smaller maximum 

absolute error in the ith state, and in that sense has outperformed the 
linear model. 

At the completion of any particular simulation session, the user has 
the option to plot the selected curves that he spooled earlier. The rou- 
tine HARDCOPY, again written in FORTRAN, reads the trajectory solutions 
from SPOOL and records at the Versatec printer /plotter , a data sheet cor- 
responding to each plot set which follows. The comparative trajectories 



17 



themselves are then plotted. In this way, permanent records of results 
are made available which are suitable for Inclusion In reports such as 
this one. 

This section has attempted to summarize in an easy-to-read format 
the functions and capabilities of the CATNAP software package. The compu- 
ting hardware situation at the University of Notre Dame has, however, at 
this time a rather significant gradient. Within the next calendar year, 
progress is expected in (1) replacing the IBM 370-168 with a newer gen- 
eration system, (2) establishing a new system in the College of Engine- 
ering, and (3) studying the possibilities for networking machines in the 
Department of Electrical Engineering. Accordingly, the specifics of 
CATNAP will have to evolve continually. Appendix B gives some details 
for students who use the package this year. While it is Illustrative, 
it can certainly be skipped on an initial reading. 


SECTION IV 
QCSEE 

Since the primary goal of the report is to apply the modeling metho- 
dology described in Section II on a practical example of a turbofan jet 
engine, it is necessary that the reader be somewhat familiar with the op- 
eration of the jet engine. Therefore, the intent of this section is to 
present an overview of the jot engine and its corresponding real-time 
digital simulator. More specifically, attention will center around NASA’ 
£uiet, Clean, JShorthaul Experimental Engine (QCSEE). 

4.1 TURBOFAN JET ENGINE 

The major purpose of the QCSEE ("Quixie") program [16,17] is to de- 
velop the technology required for quiet, clean, economical propulsion 
systems to be used in commercial short takeoff and landing (STOL) air- 
craft. As a result, QCSEE will greatly expand airline capabilities in 
densely populated metropolitan areas where low noise and atmospheric pol- 
lution requirements are critical. 

The QCSE engine, shown in Figure 4.1, will pioneer with several con- 
cepts not all currently used on existing turbofans. The engine features: 

(1) an extremely high-bypass ratio; 

(2) a gear-driven reversible-pitch fan with composite 
fan blades; 

(3) a variable geometry fan exhaust nozzle; and, 

(4) a digital electronic control system combined with 
a hydromechanical fuel control. 

The following will be an attempt to present a brief description pertinent 




Figure 4.1. QCSEE Propulsion System. 
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Co each of chese features. 

First, by doubling the ratio of cool, dense bypass air to high-temp- 
erature high-pressure air flowing through the core, QCSEE is able to re- 
duce the airflow velocity 40 7* while maintaining engine thrust, thus im- 
proving lift. Equally as important, the noise level of the jet engine is 
reduced significantly. Current estimates indicate that QCSEE will have a 
takeoff and landing "footprint" of approximately .8 square miles, whereas 
current transport aircraft have noise footprints ranging anywhere from 6 
to 80 square miles. 

Another advantage of the extremely high-bypass system is that it 
causes the engine exhaust temperature to be relatively low when compared 
to temperatures of other turbofan engines. This cooler air reduces acous- 
tic and thermal fatique in the flap systems, and therefore permits manu- 
facturers the use of conventlal airframe materials. On the other hand, 
QCSEE employs the use of exotic composite materials, whenever feasible, 
to reduce weight. Examples indicating where composite structural materials 
are used include the inlet duct and fan blades. 

The second Innovative feature of the QCSE engine is that it uses 
variable-pitch fan blades. Variable pitch will eliminate the need for 
large, complex, heavy thrust-reversal systems which are typically used on 
turbofans. The almost instantaneous response to engine thrust settings 
that variable pitch fan blades provide will give pilots improved control 
in the critical landing and takeoff portions of the flight envelope. Cur- 
rently, a 3 to 5 second time lag exists in engine acceleration from ap- 
proach thrust to full thrust. This time lag can be reduced to less than 
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one second with the QCSE engine. This is achieved by maintaining high fan 
speed, with a low thrust pitch angle. To achieve full thrust the fan pitch 
angle is merely changed accordingly. 

The next Important feature of QCSEE is its fan exhaust nozzle. The 
fan exhaust nozzle is a variable-area, four-flap design capable of area 
change from takeoff to cruise, as well as opening to a flared position to 
form an inlet in the reverse thrust mode. 

Finally, QCSEE utilizes an on-board system to control various steady- 
state and transient characteristics. The system designed to accomplish 
this task incorporates two basic control components, a modified hydromechan- 
ical fuel control and an engine-mounted, 4000 word core memory, digital 
electronic control designed specifically for QCSEE (See [18]). 

An F101 fuel pump is used in the system for supplying fuel for engine 
operation, and for operating servomechanisms in the hydromechanical control. 
A constant pressure hydraulic pump supplies fluid for operation of the ac- 
tuators which position the fan nozzle and the hydraulic motor which drives 
the variable-pitch actuators. 

An important estimation in any jet engine is the thrust. Here, the 
controller uses the ratio of compressor discharge pressure to free stream 
total pressure aS a measure of engine thrust. Incidently, the QCSE engine 
is designed to provide 17,400 pounds of installed thrust at takeoff on a 
90° fahrenheit day. 

The ultimate goal of propulsion engineers is to equip aircraft with 
fly-by-wire systems where essentially the entire flight is managed and 
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controlled by a sophisticated computer. In this arrangement, the pilot 
would merely serve as a backup. Advanced versions of QCSEE are expected 
to contain fly-by-wire systems where flight plans from one airport to an- 
other would be relayed from the ground-based center to the on-board com- 
puter. Then, during takeoff, cruise, and approach the appropriate plan 
would be initiated with the computer controlling the necessary parameters 
such as thrust, fan-blade setting, and exhaust nozzle area. 

It can be determined from the above discussion that QCSEE is a rather 
large engine made up of several bulky component systems. Consequently on 
an experimental basis, it is not very practical, let alone affordable, to 
use a real turbofan. A useful alternative, though, is to use a digital 
simulation of the engine. In this spirit, the remainder of the section 
outlines the QCSEE simulation deck which is used in the modeling exer- 
cises to follow. 

A. 2 DIGITAL SIMULATION 

With the advent of faster, more reliable computers, real-time digital 
simulation has come to play an important role in the development of ad- 
vanced aircraft and propulsion systems. Simulations are advantageous 
since, at a meager cost, they can enhance product development by elimin- 
ating the need to build expensive prototypes. Furthermore, simulation 
provides a safe means of studying various failure modes such as engine 
stall, which might otherwise be impossible to perform due to the hazar- 
dous conditions involved, Accordingly, the QCSEE simulator program [19,20] 
was established to achieve detailed real-time propulsion simulation to be 
used in manned flight simulators designed specifically for STOL aircraft. 
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For reference, Appendix A contains a flowchart of the QCSEE simulator 
as well as a list of important engine and control variables used through- 
out the simulation. 

The analytical model, developed from component performance maps of 
the various engine dynamics, incorporates multivariable functional repre- 
sentations of two or three variables. Since 40 percent of the total cal- 
culation time is consumed in function generation, a number of more effi- 
cient generating schemes are employed. Linear interpolation, polynomial 
and geometric curve fitting, and table search are among the techniques used 
to generate these functions. The eight engine dynamics chosen as state 
variables for the system include two fan speeds, four pressures, and two 
temperatures. In addition, there are four sensor states and four control 
actuator states. 

A digital controller is incorporated into the overall simulation, 
which uses absolute error signals to achieve feedback in the control loop. 
Hie error signals proceed through a segment which contains a lag compen- 
sator for stability and transient anticipation, logic for switching to 
other control modes, selectors through which limits are applied, and ele- 
ments which generate the output signals. Overides are applied to each of 
the manipulated variables under certain conditions for safety or opera- 
tional reasons. A block diagram of the QCSEE simulator control structure 
can be found in Figure 4.2. 

The computer program is written in double precision ANSI FORTRAN V, 
where great care has been taken to minimize unnecessary computations. 




Figure 4.2. QCSEE Control Structure. 
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Functions such as exponential are avoided, and whenever possible, divide 
computations are replaced by multiplications. Surprisingly, conditional 
branches are substituted for subroutine calls throughout the code. The 
program contains provisions fpr monitoring and displaying ninety engine 
and control variables, for detecting certain manfunctions, and for taking 
corrective action in the event of some critical malfunction such as fan 
drive gear failure, loss of fan speed signal, and certain control compu- 
tation faults. The integration stage of the routine utilizes a second- 
order modified Euler algorithm [21] which was found to be computationally 
faster than a standard fourth-order Runge-Kutta method. The concept of 
point averaging, rather than slope averaging, is used in this modified 
Euler algorithm, where it has been proved stable over a wide range of 
transients . 

The engine states used in the modeling discussions in the next sect- 
ion are the combustor discharge pressure (P4GS) , the core nozzle total 
pressure (P8GS) , and the rotor dynamics in the form of fan speed (NL) , 
and compressor speed (NH) . As a comparison, the fan attains a speed of 
3000 rpm, at cruise, whereas the compressor reaches a speed near 12500 
rpm. The four remaining states are the total pressures from the fan inlet 
and fan discharge, and the exit temperatures from the compressor and com- 
bustor. All three engine control inputs are employed in the modeling, 
namely, the main burner fuel flow (WFM) , the exhaust nozzle area (A18) 
and the fan pitch angle (BETAF) . To illustrate QCSEE's fuel comsumption, 
the main burner fuel flow attains a value of 4000 lb/hr at cruise, and 
a value as high as 6600 lb/hr during takeoff. 
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QCSEE has one primary input variable for manipulation, the percent 
power demand (PWRX). In reality, this power setting is a digital elec- 
tronic input supplied from the control room, or in an aircraft applica- 
tion, from the pilot or aircraft flight control computer. In the simu- 
lator, though, power demand is just a numerical percentage with idle speed 
occuring at 62. 5%. Cruise speed is considered to be about 85% PWRX, while 
full throttle occurs at 100% PWRX. A sample operation might be to run 
QCSEE into the steady state (approximately ten seconds) at idle speed, 
and then switch to a desired PWRX to perform particular transient re- 
sponse studies. Finally, the program permits modification of several 
environmental parameters, however, the simulator is normally operated at 
sea level and zero velocity. 

The complex, advanced design of the QCSE engine lends itself as a 
practical yet challenging nonlinear modeling problem worth studying. In 
the next section the QCSEE simulator is incorporated into CATNAP and used 
to formulate nonlinear models of the type mentioned thus far. 


SECTION V 


QCSEE MODELING 

A useful way to test the feasibility of the modeling technique under 
investigation is to apply the ideas presented thus far using a nonlinear 
digital simulation. Accordingly, this section summarizes a research ef- 
fort focused on developing nonlinear models of tensor type for the QCSE 
engine by using the CATNAP software discussed in Section III. Specific- 
ally, the aim of the section is two-fold. First, the problem is formu- 
lated in reference to the jet engine digital simulation. This involves 
the creation of an appropriate CATNAP loader routine which utilizes the 
QCSEE deck. Then, an example four-state, three-control nonlinear model 
of the QCSE engine is presented. The majority of this section consists 
of simulation results in the form of solution plots and tabulated error 
criteria. These test results are provided to illustrate a comparison of 
the true, linear, and nonlinear solutions, and to establish an acceptable 
region of model validity. 

5.1 FORMULATION 

Attention in this section will center around the formulation of a 
reduced order four-state, three-control analytical model. As mentioned 
in the previous section, the chosen states are the combustor discharge 
pressure (P4GS) , the core nozzle total pressure (P8GS) , the fan speed 
(NL) , and the compressor speed (NH) . Recall also that the main burner 
fuel flow (WFM) , the exhaust nozzle area (A3, 3) , and the fan pitch angle 
(BETAF) are the control inputs employed. Since the region of interest 
for identification is the operating envelope, an expansion is made about 
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a particular non-zero point, (Xq,u^), on the envelope. If the four- 
tuple x is given by the expression 

x(t) - 5(t) - x Q (5-1) 

where 


x^ " P4GS (psi) , 

X 2 " P8GS (psi) , 

X 3 - NL (rpm) , 
x 4 “ NH (rpm) , 

and the three- tuple u is taken to be 

u(t) ■ u(t) - Uq (5-2) 

where 

u. ■ WFM (lb /hr) , 

1 m 

u 2 « A18 (in 2 ) , 
u 3 ■ BETAF (deg.) , 


then from Section 2.1 the differential equation representation, with 

tensor notation applied, has the following form: 

. O' .O' 

x * L 10 X + L 01 u 


+ L 20 (x ® x) + L 1:1 (x ® u) 
+ Lq 2 (u ® u) + . . . 


(5-3) 


It should be noted that in the off-simulator examples (See [ 3 ]) the 
equilibrium point (point of expansion) was taken to be zero so that a 
multivariable equivalent of a Maclaurin series expansion was used. In 
this study, however, the point of expansion is not zero and a Taylor 
series of several variables is used. 
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The model presented in this section is of degree two; that is, 
second degree tensor terms 

(x ® x) , (x ® u) , (u ® u) 

are the highest degree terms retained. Therefore, the truncation approx- 
imation of equation (5-3) takes on the structure 

A* 1V1 Aj Aj Aj Aj 

x - [L 1q L Q1 L 2q L 1;l L q2 ] x t . (5-4) 

In this particular situation x T contains thirty-five monomials, seven 
of degree one, and twenty-eight of degree two. These are 


(x x , x 2 , x 3 , x 4 , Up u 2 , u 3 , 

*V x l x 2» x l x 3’ 

w 

x 2* X 2 X 3* X 2 X 4 ' 

X Z y X 3 X 4* x 4> KjUp x 

l u 3 , x 2 u l* X 2 U 2’ 

x 2 u 3 , 

X 3 U 1 J X 3 U 2’ 

x 3 u 3 , x 4 u , x A u 2 , x 4 u 3 , uj. 

u l u 2 , u x u 3 , ujj, 

U 2 U 3’ 

u*) T . (5-5) 


Here, the elements of (5-5) are listed according to the established or- 
dering convention. 

Appropriate engine operation for the model identification involves 
opening the loop by deactivating the controller and independently insert- 
ing the individual control inputs (See [22]). In this strategy, nonlin- 
earities of the plant exist — which might otherwise be less noticeable 
had the controller been present in the loop. This scheme is portrayed 
in Figure 5.1 where a constant desired power percentage (DPP) is applied 
and the controller is effectively "turned-of f" by equating the controller 
state derivatives with zero. In this way, an input vector of the form 

(u) ± - (u 0 ) ± + o ± cos(2tt d> i t) , i - 1,2,3 (5-6) 

can be substituted, where a and <j> are three- tuples containing ampli- 
tudes and frequencies (in hertz) respectively. An important point to 



acos (2]r<{>t) 



•pen Loop Identification Scheme. 
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note Is that, In the open loop situation, the choice of Input signals is 
critical. This is due to the fact the engine itself has certain physical 
limits, which in turn have been placed in the simulator. In reality, 
exceeding these limits could cause severe damage to the engine, an ex- 
ample of which is turbine melt down. 


Using the QCSEE digital simulation, and implementing the above ideas, 
a CATNAP loader routine (QCSELOAD) is constructed to set up the problem 
for the least squares identification. Moreover, if h time points are 


used in the sampling procedure, the matrix equation 

X “ ^10 L 01 L 20 L 11 L 02 J *T 
4 x h 4 x 35 35 x h 


(5-7) 


is formulated, with matrix dimensions as shown, for fourth order models 
keeping second degree tensor terms in the approximation. In addition, 
provisions are made available in the code which permit normalization 
of the sampled data according to the scheme outlined in Section 2.2. 


In order to perform comparative model simulations, the engine pro- 
gram must be incorporated into a CATNAP simulator routine. Consequently, 
program QIXSIM is built using QCSEE as the true system in the application 
subroutine TRUES. QIXSIM is used extensively in the next section where 
results for various simulations of the example model are presented. 

5 . 2 EXAMPLE MODEL 

To aid in the selection of input signals, a family of parametric 
plots have been constructed using QCSEE steady state data from idle 
(62. 5% PWRX) to maximum power (100% PWRX) . Figure 5.2 contains an ex- 
ample of one such steady state plot. From these plots a set of accept- 



STERDT STRTE PLOTS 
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able Input signals can be selected. Acceptable state perturbations can 
be selected In a similar fashion. Another Important feature of these 
Steady state plots Is that they suggest regions of nonlinearity. From 
Figure 3.2 it appears that In the locality of 92% power demand the engine 
Is nonlinear due to the abrupt changes in exhaust nozzle area and fan 
pitch angle, with this in mind, 92% shall be established as the operating 
point of the current example. 

The following is an overview of the particular model formulation 
using QCSELOAD. The QCSEE simulator is run, closed loop, with a 92% 
power demand for ten seconds to settle all transients. This produces 
the equilibrium value (Xq,Uq). The initial conditions thus generated 
form the point of expansion for the series truncation approximation in 
the model formulation. Within the digital simulation program the con- 
troller is disconnected by setting the control derivatives to zero. 

From the steady state plots a point (£,fl) is chosen on the engine op- 
erating line at 93% power demand. The state variables are perturbed x 
from their equilibrium values where 

x - ft - Xq . (5-8) 

Furthermore, the control variables are manipulated so that the input 
vector (eq. 5-6) is inserted, where o is given by the expression 

o - a - u Q , (5-9) 

and 

$ « (2.0, 1.0, 1.3) . 

The observed states and inputs are sampled over a six second interval; 

100 samples are evenly spaced at 0.06 seconds, and the normalized dif- 
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ference between these values and the corresponding equilibrium values, 
together with the normalized ordered tensor product terms and normalized 
state derivative values comprise the data necessary for the Identifica- 
tion. Once this is complete, the model parameters contained in matrix 

'V 

representations of the are easily computed from IDENTIFY and then 

sent to the PDP 11/60 via TRANSFER. 

Using the GENERATE segment of CATNAP, two models have been iden- 
tified: a second-degree nonlinear model, and a first-degree linear 

model. The linear model has been identified for use in comparison 
studies. Accordingly, such a nonlinear model is expected to out perform 
the linear model in a region about the point of expansion. Complete 
listings of the identified model operators for both the linear and non- 
linear models are contained in Appendix C. 

Verification of the nonlinear model involves conducting numerous 
simulations for various combinations of the parameters Xq, a, end 
comparing results against those of the true system. Frequently though, 
as is the case here, a third system is included in the comparison studies 
in the form of a linear approximation. This provides the engineer with 
a gauge to determine how well a new nonlinear model performs as compared 

to a traditional representation. Utilizing the CATNAP capabilities of 

» 

QIXSIM and HARDCOPY, several validation studies have been completed to 
date, all yielding satisfactory results. 

Four such tests will be presented here, the first of which con- 
sists of twenty-five simulations along the line of engine operation with 
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different choices of x^, and a, and one triplet ($.,,» $ 2 » ^ 3 )* Re “ 
suits of this test are tabulated in Table 5.1, where 4 ■ (2.1, 1.0, 1.5) 
for all simulations. By observing the table, one can see that the error 
criterion (e^) discussed in Section 3.3 is negative for each state (x^) 
throughout all Individual simulations. This indicates that the model has 
outperformed the linear approximation in a region about the expansion 
point. To illustrate the tracking ability of the model, consider a rep- 
resentative comparison plot of these solutions against the true solution, 
Take as an example the simulation with conditions from Table 5.1, line 18 
which is depicted in Figures 5.3 and 5,4. Figure 5.3 is presented merely 
to give the reader an opportunity to see the format of a CATNAP generated 
data sheet, whereas the actual trajectory comparisons are presented in 
Figure 5.4. Clearly, the model solution tracks the true solution well 
in all four states throughout the time interval. 

The second test represents on the line studies of model performance 
for a variety of input frequencies and zero initial conditions. Table 
5.2 shows the results of these error analyses. As in the previous test, 
the are negative for all simulations. A comparison plot of one of 

these tests. Table 5.2 line 8 , is given in Figure 5.5a for variable P4GS 
and in Figure 5.5b for variable NL, The next figure shows the ability 
of the model to track a step response. Figures 5.6a and 5.6b offer 
comparison plots with conditions from Table 5.2 line 7 for variables 
P4GS and P 8 GS respectively. In both Figures 5.5 and 5.6 the model 
simulation follows the true solution well while the linear approximation 


is less accurate. 


Table 5.1: Comparison Studies on the Operating Line; ♦ * (2.1, 1.0, 1.5) 
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Figure 5.3. Data Sheet for Figure 5.4. 


COMPARATIVE SOLUTION PLOTS: STATE «* 1 

RUE A - LINEAR MODEL * - NONLINEAR MODEL 







COMPARATIVE SOLUTION PLOTS: STATE w 2 

TRUE A - LINEAR MODEL * - NONLINEAR MODEL 
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COMPARATIVE SOLUTION PLOTS: STATE * l 

RUE A - LINEAR MODEL * - NONLINEAR MODEL 





COMPARATIVE SOLUTION PLOTS: STATE # 3 

TRUE A - LINEAR MODEL * - NONLINEAR MODEL 



COMPARATIVE SOLUTION PLOTS: STATE » 1 

RUE A - LINEAR MODEL * - NONLINEAR MODEL 



TIME ISECI 

Figure 5.6a. Table 5.2, Line 





COMPARATIVE SOLUTION PLOTS: STATE » 2 

RUE A - LINEAR MODEL * - NONLINEAR MODEL 



TIME ISEC1 

Figure 5.6b. Table 5.2, Line 
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The third test considers several possible combinations of x^, and 
and fixes a ■ (55.78, -14.03, -0.198). These results are recorded 
in Table 5.3. Again, the error analyses indicate that the model outper- 
forms the linear approximation for the various test conditions. Figure 
5.7 offers an illustration which uses the conditions from line 20 of 
Table 5.3* 

The final test represents analyses for various amplitudes and fre- 
quencies less initial conditions. The unique feature of this test is 
that the excursions are made away from the typical engine line of opera- 
tion. This is done because in practice such a model may be required to 
perform within a tubular region, with the operating line at its center, 
instead of just along the line itself. The error analyses associated 
with these simulations are contained in Table 5.4. The final plots to 
be presented here are for a simulation with conditions from line 6 of 
Table 5.4. Figure 5.8a offers the comparative plot for variable P8GS 
while Figure 5.8b illustrates variable NH. 


Comparison Studies on the Operating Line; a = (55.78, -14.03, -0.198) 
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COMPARATIVE SOLUTION PLOTS: STATE * 1 

RUE A - LINEAR MODEL * - NONLINEAR MODEL 







COMPARATIVE SOLUTION PLOTS: STATE « 2 

RllE A - LINEAR MODEL * - NONLINEAR MODEL 



TIME (SEC) 

Figure 5.7b. Table 5.3, Line 20 






TIME 

Tab] 




TIME (SEC1 

Figure 5.7d. Table 5.3, Line 20 
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SECTION VI 


CONCLUSIONS 

Based upon the results included in this report it can be said that 
the tensor approach, used to model the QCSEE digital simulation, has 
tested cut well. This is true in the sense that (1) initial conditions 
and input signals have been chosen so that the nonlinear model outper- 
forms the linear approximation locally and also expands the region of 
desirable model behavior; (2) suitable nonlinear model performances 
have been observed over a wide variety of initial conditions, input 
amplitudes, and input frequencies; and (3) acceptable nonlinear model 
trajectories have been viewed in a qualitative , fashion. Besides the QCSEE 
modeling effort, important milestones have been reached, as a result of 
this study, in the area of specialized computer software development. 

With the inception of the CATNAP package, tensor nonlinear modeling 
software capability has been substantially advanced. 

A potentially useful idea which has stemmed from this work is the 
concept of nonlinear model identification using a "layered" tensor ap- 
proach. Because of the possibilities, it is perhaps worth introducing 
this approach here. 

Basically, the idea requires that the model structure (equation 
2-3) be 'jplit into several smaller identification problems. Initially, 
a first-degree approximation is formulated, using the original techni- 
que, and identified. This is called the first layer. Next, a second- 
degree model is formulated; however, this time the response from the 
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first-degree model is constructed and subtracted from the derivative es- 
timate matrix of the second-degree formulation. The second-degree form- 
ulation uses the identified first-layer parameters and the appropriate 
first degree monomials. This results in a second-layer formulation 
which contains only second degree tensor terms. Once identified, these 
parameters (second-layer) are combined with the first-layer parameters 
to form a model structure containing both first and second degree tensor 
terms. It should be noted that the magnitude of the initial condition 
and input perturbations used for excitation usually increase as the 
layer number does. This is done to assure that the derivative estimates 
maintain sufficient rank. 

To add third degree tensor terms, the second-degree model response 
is constructed, and subtracted from the derivative estimate matrix of 
the third-degree formulation. A third-layer formulation containing 
only third-degree tensor terms is the result. When the third layer 
parameters are identified they are combined with the first and second 
layer ones to form a model structure. Clearly, any number of layers 
could be added using this scheme. 

Intuitively, one would expect this approach to reduce the size of 
the nonlinear model identification problem significantly, as well as 
reduce the chance of obtaining dependent data. In addition, this ap- 
proach should enlarge the region of acceptable model behavior. In this 
regard, it is believed that further studies are warranted. In particu- 
lar, a thorough investigation of the layered tensor approach is recom- 
mended. Some preliminary efforts have already begun in this area. 


APPENDIX A 

QCSEE SIMULATION FLOWCHART AND VARIABLE LIST 
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FLOWCHART FOR QCSEE DIGITAL SIMULATION 
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ALT 

MNO 

PO 

TO 

VO 

FN 

FGROSS 

FRAM 

FG8 

FGH 

FG18 

Pll 

PS11 

XM11 

P12Q11 

P12GDT 

P12GS 

P13GDT 

P13GS 

P13Q12 

P18 

P25 

P3 


VARIABLE LIST FOR OCSEE DIGITAL SIMULATION 

- Altitude - ft. 

- Mach Number ~ dim. 

- Ambient Pressure - psia 

- Ambient Temperature - *R 

- Frees tream Velocity 

- Engine Net Thrust - lb^ 

- Engine Gross Thrust - lb f 

- Engine Ram Thrust ~ lb£ 

- Core Engine Gross Thrust ~ lb^ 

- Fan Reverse Gross Thrust ~ lb^ 

- Fan Forward Gross Thrust ~ Ibj 

- Engine Inlet Total Pressure - psia 

- Engine Inlet Static Pressure - psia 

- Inlet Mach Number dim, 

- Inlet Total Pressure Ratio, P i2^ P ll ~ 

- Fan Inlet Duct Total Pressure Iteration Potential ~ psia/sec 

- Fan Inlet Duct Total Pressure Updated Iteration Value ~ psia 

- Fan Discharge Total Pressure Iteration Potentis.A ~ psia/sec 

- Fan Discharge Total Pressure Updated Iteration /alue ~ psia 

- Fan Pressure Ratio - P^/P^ 

- Fan Nozzle Total Pressure - psia 

- Compressor Inlet Total Pressure - psia 

- Compressor Discharge Total - psia 

- Compressor Pressure Ratio - P^/^ 5 




P3Q25 
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PS3C 

EPRSN 

P4GDT 

P4GS 

P42 

P8GDT 

P8GS 

T12 

T13 

T25 

T3 

WQW12 

W12R0 

W12 

W18 

W25 

W3 

W4 

W8 

T4 

T41C 

T42 

T8 

HPF 

HPLT 


- Compressor Discharge Pressure at Sensor Location * psia 

- Sensed Engine Pressure Ratio ~ dim. 

~ Combustor Discharge Pressure Iteration Potential ~ psia/sec 

- Combustor Discharge Pressure Updated Iteration Value ~ psia 

- Low Pressure Turbine Discharge Total Pressure - psia 

- Core Nozzle Total Pressure Iteration Potential ~ psia/sec 

- Core Nozzle Total Pressure Updated Iteration Value ~ psia 

- Fan Inlet Temperature - *R 

- Fan Discharge Temperature - °R 

- Compressor Inlet Temperature - *R 

- Compressor Discharge Temperature - °R 

l 

- Fan Flow Function Ratio ~ dim. 

- Forward Fan Map Flow Function ~ dim, 

- Fan Airflow - lb /sec 

HI 

- Fan Nozzle Airflow - lb /sec 

m 

- Compressor Inlet Airflow - lb /sec 

IB 

- Compressor Discharge Airflow - lb m /sec 

- Combustor Airflow - Ih/mc 

tn 

- Core Nozzle Airflow - lb /sec 

m 

- Combustor Discharge Temperature - °R 

- Calculated Turbine Inlet Temperature - °R 

- High Pressure Turbine Discharge Temperature - °R 

- Low Pressure Turbine Discharge Temperature - °R 

- Fan Horsepower ~ hp. 

- Low Pressure Turbine Horsepower ~ hp. 
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NLDT 

NL 

HPC 

HPHT 

NHDT 

NH 

PCNLR 

PCNHR 

PCNHRD 

PCEPRD 

PWRX 

ICTRMO 

IB1F 

NLK 

NLMXR 

XMV 

WFM 

X18 

A18 

EPRD 

DFSEL 

EPWR 

ETWR 

ENWRH 

ENWRL 


- Fan Speed Derivative - rpm/sec 

- Fan Rotor Speed - rpra 

- High Pressure Compressor Horsepower ~ hp. 

High Pressure Turbine Horsepower ~ hp. 

- Compressor Speed Derivative - rpw/sec 

- Compressor Speed - rpra 

- Percent Corrected Fan Speed 

- Percent Corrected Core Speed 

- Percent Corrected Compressor Speed Demand 

- Percent Demanded Engine Pressure Ratio 

- Power Demand - Percent 

- Control Mode Logic Variable, 0 - Forward, 1 ■ Reverse 

- Inadvertent Reverse and Out-Reverse Interlock Logic Signal 

- Scheduled Fan Speed Demand ~ rpm 

- Fan Speed Selector Value ~ rpm 

- Fuel Metering Value Position - in. 

- Main Burner Fuel Flow - lb /hr 

m 

- Fan Nozzle Area Actuator Position - in. 

2 

- Fan Nozzle Area - in 

- Demanded Engine Pressure Ratio ~ dim. 

- Digital Fuel Control Selector ~ in/sec 

- Digital Engine Pressure Ratio Control Error ~ In/sec 
-Digital Maximum Turbine Inlet Temperature Control Error ~ in/sec 

- Digital Core Speed Control Error ~ in/sec 

- Digital Fan Speed Control Error ~ in/sec 
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ENWRI 

QMV 

QEOR 

QNH 

QDEC 

QACC 

BETAF 

DXSEL 

EX18 

EX18MX 

EX18MN 

EX18RV 

DBSEL 

EBNBR 

EBMR 

EBREV 

THETAI 

BFW 


- Digital Core Idle Speed Control Error ~ in/sec 

- Hydromechanical Power Piston Flow ~ in / sec 

3 

- Equivalent Electronic Override Power Piston Flow ~ in /sec 

3 

- Hydromechanical Core Speed Control Power Piston Flow ~ in /sec 


- Hydroraechanical 
~ in 3 /sec 


Deceleration Limit Control Power Piston Flow 


Hydroraechanical 
~ in 3 /sec 


Acceleration Limit Control Power Piston Flow 


- Fan Pitch Angle - Degree 

- Digital Fan Nozzle Position Control Selector ~ in/sec 

- Digital Inlet Mach Number Control Error ~ in/sec 

- Digital Maximum Fan Nozzle Position Control Error ~ in/sec 

- Digital Minimum Fan Nozzle Position Control Error ~ in/sec 

- Digital Fan Nozzle Reverse Position Control Error ~ in/sec 

- Digital Fan Pitch Control Selector ~ in/sec 

- Digital Fan Pitch Angle Floor Schedule Control Error ~ V/V/sec 

- Digital Maximum Fan Pitch Angle Control Error ~ V/V/sec 

- Digital Reverse Fan Pitch Angle Control Error ~ V/V/sec 

- Fan Pitch Mechanism Drive Motor Position ~ deg. 

- Fan Pitch Position Sensor Signal ~ V/V 
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APPENDIX B 


CATNAP USER’S GUIDE 

This appendix takes the user step-by-step through the operation of 
CATNAP to identify and simulate a nonlinear model. The process is divided 
into three main sections and covers everything from logging-on to logging- 
off. 

The first section is devoted to identifying a nonlinear model and 
transfering it to the Radiation Laboratory using the GENERATE segment of 
CATNAP which is resident on the IBM 370-168 computer system. The process 
of retrieving the nonlinear model from the Radiation Laboratory and enter- 
ing it into the model library on the electrical engineering department 
PDP 11/60 computer is the concern of the second section. The third and 
final section uses the CATNAP segment SIMULATE, which is on the PDP 11/60, 
and leads the user through a sample simulation session. 

B.l MODEL IDENTIFICATION PROCEDURE 

This section takes the user through the process of identifying a non- 
linear model. This step-by-step process begins with logging-on the IBM 
370-168, checking the user library, entering the GENERATE supervisory, 
executing a desired loader routine, executing the SPEAKEASY program 
IDENTIFY, sending the model to the Radiation Laboratory and finally 
logging-off the system. 

To log on, the user must initiate a connection between the terminal 
and the computer, and then identify himself to the system as an authorized 
user by entering a valid TSO identification number. The discussion to 
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follow will be directed to users of DECwriter terminals only, since the 
procedure for making the connection is terminal dependent. 

To log on at a DECwriter, the user must turn the power switch on 
the keyboard to ON, and check to make sure that the two switches on the 
back of the telephone coupler are in the power on and half duplex mode. 

On the keyboard, the CAPS LOCK key should be in the "down" position. In 
addition, the panel to the left of the keyboard should have all keys in the 
"up" position except for the HDX (half-duplex) key and the 300 BAUD key, 
which should be down. Lift the telephone receiver, dial the remote ac- 
cess telephone number and listen for a tone. Upon hearing the tone, place 
the receiver in the coupler with the cord toward the back of the coupler. 
When the carrier light on the front of the coupler becomes illuminated 
strike the capital letter "X" a couplfe of times and then hit the carriage 
return key (referred to from now on as CR) on the right side of the key- 
board. (Note: Throughout this appendix, information printed by the com- 

puter will be underlined; all other information is user input.) When the 
proper connection is made, the terminal will print a log on message. 

IKJ54012A ENTER LOGON- 

Respond by typing 

LOGON XXXXXX SIZE (1000) (CR) 

where XXXXXX is the user’s TSO ID and the SIZE keyword requests addi- 
tional memory for the session. If the system accepts the ID it will 
print an acknowledgement. 

' XXXXXX LOGON IN PROGRESS 


The computer will then print any broadcast messages and the word READY. 
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READY 

It is usually a good idea to list the user library to see what is 
stored in it. The user should have on hand the following files in order 
to identify models: CLIST, CNTL, SPEAKEZ.DATA and MODEL. DATA. Also, 

the user should have on file a FORTRAN application loader routine for 
the particular problem being studied. The SPACELST command is used to 
list the files stored. 

SPACELST (CR) 

The computer responds by listing the names of all files stored on 
disk under the user’s TSO ID. 


OTHER FILES 

■ 

CLIST 

CNTL 

MODEL. DATA 
QCSELOAD . DATA 
SPEAKEZ.DATA 

READY 

Seeing that the required files are present, the user can now access the 
GENERATE program. Incidently, QCSELOAD. FORT is the loader routine used 
to formulate nonlinear models of the QCSEE simulator. 

To enter the supervisory the user must execute the proper command 
list program. 

EX (GENERATE) (CR) 

The program will begin to run and respond by printing an introduction. 
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WELCOME TO THE NONLINEAR MODEL GENERATION S EGMENT OF THE COMPUTER 
AIDED TENSOR NONLINEAR MODELING APPLICATIONS PACKAGE (CATNAP) t 

ENTER NAME OF DESIRED LOADER ROUTINE ; 

The user can then enter the name of the loader routine to be executed. 
As an example, say the user wishes to run QCSELOAD . FORT . The user would 
then merely enter QCSELOAD as follows; 

ENTER NAME OF DESIRED LOADER ROUTINE; QCSELOAD (CR) 

Control would then pass from the supervisor to QCSELOAD and execution 
would begin. (Note; Since the purpose of this appendix is to outline the 
CATNAP package in general and not focus on a particular application rou- 
tine, QCSELOAD will not be addressed here; however, Appendix C contains 
listings of sample model formulations and identifications which utilize 
QCSELOAD.) 

At the termination of the loader routine control is passed to the 
SPEAKEASY processor and the message 

TSO SPEAKEASY III PI 

is printed at the terminal followed by the prompt symbol (;__). The user 
then types the string 

;_ SIZE- 500 ; GET - IDENTIFY : IDENTIFY : QUIT (CR) 
to begin program IDENTIFY which prints the model coefficients at the 
terminal and in file MODEL. DATA. Once the model coefficients are printed 
control is returned to the supervisor where a couple of questions are 
asked . 

DO YOU WISH TO SEND THIS MODEL TO THE RADIATION LABORATORY? fY/N] ; 


A "Y" response to this question would initiate execution of program 
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TRANSFER which would send a copy of MODEL. DATA Co the PDP 11/55 computer 
located in the Radiation Laboratory. An "N” response would merely cause 
the supervisor to print the next question. 

DO YOU WISH TO IDENTIFY ANOTHER MODEL? [Y/N] ; 

A ,f Y M response here would cause the program to branch back to the question 
ENTER NAME OF DESIRED LOADER ROUTINE; 
whereas an M N H response would cause the supervisor to exit. 

READY 

Once in this mode, the user is ready to disengage connection with the 
IBM 370-168. This is achieved by entering the NDLOGOFF command. 

NDLOGOFF (CR) 

The computer will respond with 

XXXXXX LOGGED OFF TSO 

and the carrier light will turn off. The user should then hang up the 
telephone and turn the power switches on the coupler and DECwriter off. 

B.2 NONLINEAR MODEL RETRIEVAL 

The procedure to retrieve nonlinear models from the Radiation Labor- 
atory is the topic presented in this section. 

Before beginning, the user must select the correct RK05 disk pack 
from the disk cabinet in the EE department PDP 11/60 laboratory. This 
disk pack is labeled CNTLl. It should be placed in a plastic bag (also 
located in the cabinet) in order to prevent the disk from collecting 
dust or getting wet while in transit. The disk is now ready to be carried 
to the Radiation Laboratory. 
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Upon entering the FDP 11/55 computer laboratory In the basement of 
the Radiation Laboratory the user should remove the disk pack from the 
bag and place it in one of the two RK05 disk drives available (either 
DK0: or DK1:). Press the RUN/LOAD switch in front and note that the RDY 
(ready) light should turn on. Also, if the WTPROT (write protect) lamp 
is lit, the WTPROT switch should be depressed so that the light goes out. 
Now the disk is loaded and the user is free to log on. 

When one of the available terminals is turned on the computer sends 
a prompt (>) and waits for the user to supply a command. To log on the 
user types HELLO. 

> HELLO (CR) 

Next the computer prompts for the user’s account number (UXC) or 
last name. For this discussion the correct UIC is [45,1]. 

>ACCOUNT OR NAME; [45.1] (CR) 

To complete the log on procedure the user's password is entered. Notice 
that in order to maintain privacy the password itself is not echoed back 
to the terminal. 

>PAS SWORD; (password) (CR) 
not echoed 

Following any broadcast messages the computer will respond with the prompt 
symbol. 

It is now necessary to allocate and mount the disk for use. To 

achieve this the user would enter 

>ALL DKn: (CR) 

>MOU DKn: CNTL1 (CR) 
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where DKn; is either DK0: or DK1: depending upon which RK05 drive is 
being used. From this point onward let us assume that DK1; is being 
used. 

Once this is complete the user has access to his disk and is ready 
to begin the copy procedure. Prior to performing the copy, the user 
must list his directory to find out the name of the model file, This 
involves using the PIP command. 

> PIP /LI <CR) 

The computer will respond by listing all the names of files stored on 
the default disk under the user’s UIC. 

OTHER FILES 

9 

• 

RADLAB.437 

777.437 

The two files of interest are RADLAB and 777. RADLAB contains the JCL 
code used to send the model over and is no longer of any use, whereas 
777 contains the nonlinear model coefficients and is the file we wish 
to retrieve. 

It is important to mention here that the ,437 which is appended 
to the file name indicates that RADLAB and 777 were contained in the 
437th job sent over this particular remote job entry port on a given 
day. This is useful to know. For example, say three models were con- 
secutively sent to the Radiation Laboratory during an identification 
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session. Then, three RADLAB files and three 777 files would be present 
on the PDP 11/55. Distinguishing then from each other might be difficult 
without the job number; however, with the appended number the order of 
arrival becomes obvious. 

Once the name of the desired model file is known, the copy can be 

initiated via a set of PIP commands. An example Illustrating this for 

file 777.437 might be as follows: 

_>PIP DK1: [45,1]*.*;*/DE (CR) 

>PIP DK1: [45, ljQ92S4MA.DAT; 1*77 7. 437 (CR) 

>PIP RADLAB. 437;1, 777. 437;1/DE (CR) 

The first line deletes any files on DK1:[45,1] that are present from a 

previous transfer session. Line two performs the actual copy, and the 

/ 

last line deletes the necessary files on the default disk. It should be 
noted that the name Q92S4MA is chosen to be unique of other model names 
in the library, and furthermore is the name that will be requested later 
in the SIMULATE segment. (Comment: Q92S4MA is the name of the QCSEE 
model discussed in Chapter V.) 

When the user has completed these steps he logs off the computer by 
typing BYE: 

>BYE (CR) 

The terminal is now ready to be turned off. 

Moving over to the disk drive, the user should depress the RUN/LOAD 
switch. Upon doing this the RDY light will turn out. The user MUST 
wait approximately two minutes until the LOAD light comes on. Only then 
can he remove the disk. from the drive and put it back into the 
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plastic bag. (Warning: Removing the disk prematurely will sever ly damage 

the disk drive and disk pack.) 

Upon returning to the PDP 11/60 laboratory in the engineering build- 
ing the disk should be loaded in the drive according to the same proce- 
dure as before. The user now logs on one of the Tektronix terminals. 

^HELLO (CR) 

This time however, the UIC is different but the password is the same. 

> ACCOUNT OR NAME: [17,1] (CR) 

>PA3SW0RD: (password) (CR) 

""not echoed 

When the user is logged on, a short indirect coramand file (ICF) is auto- 
matically executed. Its purpose is to allocate and mount the disk as 
well as program the various special function keys on the Tektronix termi- 
nal. After it is finished the prompt symbol (>) will appear. 

To complete the model retrieval process (for the example presented), 
the user types the command 

>PIP DK2: [17, 7]Q92S4MA.DAT;l-DK2s [45, 1JQ92S4MA.DAT (CR) 

which enters the file Q92S4MA into the model library located in directory 

DK2: [17,7] . 

The user is now free to log off or proceed to the next section * 

B, 3 MODEL SIMULATION PROCEDURE 

This final section takes the user step-by-step through the process 
of simulating a nonlinear model on the PDP 11/60. If the user is not 
already logged-on at a Tektronix terminal he should do so by entering 
the following commands. 


77 


>HELLO (CR) 

^ACCOUNT OR NAME: [17,1] (CR) 

>PASSWORD i (password) (CR) 
not echoed 

When the log on Is complete, that Includes the automatic execution of the 
LOGIN ICF, and the user sees the prompt symbol (>) he should strike and PT 
(point) key to the right of the keyboard followed by a CR. ^SIMULATE will 
appear on the screen and execution of the CATNAP segment SIMULATE will be- 
gin. Figure B.l contains a diagram Indicating the prompt sequences avail- 
able in the supervisory program. To illustrate the use of SIMULATE, con- 
sider the following sample simulation session. 

> Q SIMULATE (CR) 

The screen is cleared and an introduction is displayed. 

>1 WELCOME TO THE NOTRE DAME DEPARTMENT OF ELECTRICAL ENGINEERING 
>± COMPUTER-AIDED TENSOR NONLINEAR MODELING APPLICATIONS PACKAGE 1 
>; 

>* £2. X2£ WISH TO EXECUTE AN EXISTING SIMULATION ROUTINE? [Y/N] : Y (CR) 

>* ENTER NAME OF DESIRED SIMULATION ROUTINE : QIXSIM (CR) 

(Note: QIXSIM is the QCSE engine simulation routine which incorporates 

the QCSEE digital deck in subroutine TRUES.) The screen is cleared again 

and program QIXSIM starts executing. 

WELCOME TO THE DYNAMICAL SYSTEM SIMULATION ROUTINE OF THE COMPUTER 
AIDED TENSOR NONLINEAR MODELING APPLICATIONS PACKAGE (CATNAP) ! 

ENTER £ OF STATES & £ OF CONTROLS TO BE USED THROUGHOUT THIS SESSION: 

The user responds by typing: 

4 3 (CR) 

The program then asks a series of questions which, of course, the user 


answers . 
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< 


Figure B.l. Prompt sequence diagram for the SIMULATE supervisory program. 
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SHOULD THIS SESSION INCLUDE THE USE OF A TRUE SOLUTION ROUTINE? fY/Nl : 

Y (CR) 

SHOULD THIS SESSION INCLUDE THE USE OF A LINEAR MODEL? fY/Nl: 

Y (CR) 

ENTER NAME OF DESIRED LINEAR MODEL FILE: 

Q92L4MA (CR) (QCSEE 4 state linear model at 92% power demand) 

USING LINEAR MODE L FILE: DK2; [17 ,7 1Q92L4MA.DAT 

SHOULD THIS SESSION" INCLUDE THE USE OF A NONLINEAR MODEL? fY/N): 

Y (CR) 

ENTER NAME OF DESIRED NONLINEAR MODEL FILE: 

Q92S4MA (CR) (QCSEE 4 state nonlinear model at 92% power demand) 

USING NONLINEAR MODEL FILE: DK2: [17 ,7 1Q92S4MA.DAT 

The above set of questions is used to establish the session configuration 
as well as isolate the necessary model files. 


Following this, the screen is cleared once more and a problem sum- 
mary is displayed. 

****************PRC)BLEM SUMMARY**************** 
CONFIGURATION: TRUE, LINEAR & NONLINEAR 
NUMBER OF STATES: 4 
NUMBER OF CONTROLS: 3 
LENGTH OF TENSOR TERM VECTOR: 35 
DEGREE OF APPROXIMATION: 2 
*********************************************** 


Next the user is prompted for various simulation parameters. 

ENTER INTEGRATION STEPSIZE, & UPPER TIME LIMIT OF INTEGRATION: 

(answers) (CR) 

ENTER THE 4 INITIAL CONDITIONS : 

(answers) (CR) 

ENTER AMPLITUDE, & FREQUENCY . . .£F FORCING FUNCTION FOR CONTROL // 1: 
(answers) (CR) 

ENTER AMPLITUDE, & FREQUENCY . . .OF FORCING FU NCTION FOR CONTROL //2: 
(answers) (CR) 

ENTER AMPLITUDE , & FREQUENCY . . .OF FORCING FUNCTION FOR CONTROL // 3: 
(answers) (CR) 

ENTER DESIRED POWER SETTING: 

92.0 (CR) 

Once these questions are answered the simulator begines to integrate the 
three systems. It should be mentioned that the last question regarding 
power setting is specific to this problem and is contained in the appli- 
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cation subroutine TRUES Instead of the mainline as are the other questions. 

After the solutions are sampled and the error criteria are computed 
and displayed the user Is asked: 

DO YOU WISH TO HAVE THE SOLUTIONS PRINTED AT THE VERSATEC? [Y/N] : 

For most simulations, an "N" response here is appropriate since the raw 
data is not really very useful. A more useful output form is the compara- 
tive trajectory plots which are displayed at the Tektronix terminal when 
the user answers "Y" to the next question. 

DO YOU WISH TO HAVE THE SOLUTIONS PLOTTED AT THE TEKTRONIX 4025? TY/Nl : 

If the user wants to make hardcopy versions of these plots he may wish to 
answer "Y" to the following question. 

DO YOU WISH TO HAVE THE SOLUTIONS FILE WRITTEN FOR HARDCOPY PLOTTING? TY/Nl : 
An "N" response would cause the program just to skip to the next question. 
DO YOU WISH TO RERUN SIMULATE WITH THE CURRENT CONFIGURATION? [Y/N] : 

A "Y" response to this prompt would cause the program to jump back to 
ENTER INTEGRATION STEPSIZE, & UPPER TIME LIMIT OF INTEGRATION: 
whereas an "N" response would cause the program to end; and control 
would be returned to the supervisor, 

When the supervisor receives control it initiates yet another prompt. 
>* DO YOU WISH TO SPOOL THE HARDCOPY PLOT FILE? [Y/N]: 

If the user answered H Y" to the earlier question regarding hardcopy plot- 
ting, then he should answer "Y" here. An "N" response at this time term- 
inates CATNAP. When a "Y" answer is detected, the program HARDCOPY is 
executed which prints corresponding data sheets along with the desired 
plots at the Versatec. 
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WELCOME TO THE HARDCOPY PLOTTING ROUTINE OF THE COMPUTER-AIDED 
TENSOR NONLINEAR MODELING APPLICATIONS PACKAGE (CATNAP) 

SHOULD DATA SHEETS BE SEPARATED BY A BLANK PAGE?[Y/N]: 

A "Y" response to this question is only necessary when report quality 

data sheets are required; otherwise, "N" should he entered to reduce 

paper consumption. 

When more than one curve is present on a plot the user is given the 
choice of having single vertical axis scaling or multiple vertical axis 
scaling. 

SHOULD THE PLOTS HAVE SEPARATE VERTICAL SCALES FOR EACH CURVE? fY/N] : 

Here, an "N" response is generally preferred since it Implements the 
more desirable single vertical axis scaling. 

Finally, the data sheets and plots are printed at the Versatec. At 
the completion of program HARDCOPY, CATNAP terminates and the user is free 
to log off. 

> BYE (CR) 
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